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Abstract 

We propose a new unified evolution equation for parton distribution func- 
tions appropriate for both large and small Bjorken variables x, which is an 
improved version of the Ciafaloni-Catani-Fiorani-Marchesini equation. In 
this new equation the cancellation of soft divergences between virtual and 
real gluon emissions is explicit without introducing infrared cutoffs and next- 
to-leading contributions to the Sudakov resummation can be included sys- 
tematically. It is shown that the new equation reduces to the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi equation at large x, to the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation at small x, and to the modified BFKL 
equations with a dependence on momentum transfer Q and with unitarity, if 
hard virtual gluon contributions and correction to strong rapidity ordering 
are retained, respectively. 
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1. Introduction 



It is known that the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation |l| sums large logarithmic corrections In Q to a parton distribution 
function at large x, with Q being momentum transfer and x the Bjorken vari- 
able, and that the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation B sums 
ln(l/x) at small x. The Ciafaloni-Catani-Fiorani-Marchesini (CCFM) equa- 
tion which is appropriate for both large and small x, has been proposed to 
unify the above two equations. The conventional derivation of the evolution 
equations involves all-order summation of ladder diagrams with rung gluons 
obeying specific kinematic orderings. For the DGLAP, BFKL, and CCFM 
equations, they are transverse momentum ordering, rapidity ordering, and 
angular ordering, respectively. 

Recently, we have applied the Collins-Soper (CS) resummation technique 
[|j] to all-order summation of various large logarithms |5| . This technique was 
developed originally to organize double logarithms In Q. We have demon- 
strated that it is also applicable to single-logarithm cases, i.e., the derivation 
of the evolution equations mentioned above ||. In this approach all the 
evolution kernels are obtained from the same one-loop diagrams, which are 
evaluated under appropriate soft approximations that correspond to specific 
kinematic orderings. Therefore, our method is simpler, and provides a unified 
viewpoint to the evolution equations. It is easy to modify the approximation 
employed in the evaluation of the BFKL kernel, such that a Q dependence 
can be introduced into the BFKL equation 0, which is important for the 
explanation of experimental data, and that predictions from the modified 
BFKL equation satisfy the unitarity bound J7|. 

As reproducing the CCFM equation in the framework of resummation, 
we notice that real and virtual gluon emissions must be treated in a differ- 
ent way: the former are computed to lowest order, generating a splitting 
function, while the latter are resummed to all orders, giving a Sudakov form 
factor. This observation is consistent with that in the conventional deriva- 
tion || . However, infrared divergences must cancel between real and virtual 
gluon emissions order by order. Hence, in the conventional derivation in- 
frared cutoffs are necessary for regularizing the individual divergences of real 
and virtual corrections. These cutoffs may cause an ambiguity, when one 
intends to include next-to-leading logarithmic summation into the evolution 
equation. 
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In Ref. || we have proposed a new unified evolution equation, which 
results from one of the applications of the CS formalism to logarithmic sum- 
mations. In this paper we shall present a detailed derivation of the new 
equation, explain the ingredients that are added compared to the CCFM 
equation, and discuss predictions for the gluon distribution function. As 
stated above, infrared divergences of real and virtual gluon emissions should 
cancel exactly without introducing artificial cutoffs, and next-to-leading con- 
tributions to the Sudakov resummation can be included systematically. To 
achieve these goals, real gluon correction is separated into two pieces. The 
piece containing an infrared pole is combined with the corresponding virtual 
correction, both of which are resummed into a Sudakov form factor. The 
other piece, being infrared finite, is evaluated to lowest order, giving a split- 
ting function. With such an arrangement, both the Sudakov form factor and 
the splitting function are rendered infrared finite. Using the CS formalism, 
the running of coupling constants appearing in the Sudakov form factor and 
in the splitting function can be taken into account. 

It will be shown that hard virtual gluon contributions introduce a de- 
pendence on the kinematic variable Q to the rise of the gluon distribution 
function at small x. This Q dependence is essential for the explanation of 
experimental data as shown in ||. In this respect, the new equation, with 
the Q dependence arising from the In Q summation, is similar to the CCFM 
equation. It has been observed in that the assumption of strong rapidity 
ordering overestimates real gluon contributions. This is the reason predic- 
tions from the BFKL equation violates unitarity constraints. In the deriva- 
tion of the new equation we shall not assume rapidity ordering, and obtain a 
destructive correction to the BFKL rise. With this correction, the increase 
of the gluon distribution function becomes mild, and satisfies the unitarity 
bound. The underlying mechanism of the above features are basically the 
same as that of the modified BFKL equations with a Q dependence and with 
unitarity || [7|. We shall argue that predictions from the CCFM equation 
violate unitarity constraints. 

The CCFM equation is reproduced using the CS resummation technique 
in Sec. 2. The new unified evolution equation is derived in Sec. 3. We demon- 
strate in Sec. 4 that the new equation reduces to the DGLAP equation, to 
the BFKL equation, and to the modified BFKL equations with a Q depen- 
dence and with unitarity under different approximations. Section 5 is the 
conclusion. 
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2. CCFM Equation 



We study the unintegrated gluon distribution function defined in the axial 
gauge n ■ A = 0, n being a gauge vector: 

x^Mtfto/-, yT)F^(0)\p,a) , (1) 

where \p, a) denotes the incoming hadron with light-like momentum p^ = 
p + 5 tl+ and spin a, and F£ is the field tensor. An average over color is 
understood. F(x,kT,p + ) describes the probability of a gluon carrying a 
longitudinal momentum fraction x and transverse momenta fey. 

For = lying on the light cone, Eq. (|1]) is gauge invariant. This 
can be understood by considering a variation of n on the light cone, under 
which Eq. ([I]) is still defined for A + = and remains the same. To derive an 
evolution equation of a parton distribution function using the CS formalism, 
we must allow n to vary arbitrarily away from the light cone (n 2 ^ 0), 
and a parton distribution function becomes gauge dependent. However, it 
will be found that an evolution kernel turns out to be independent of n. 
This is natural, since it has been proved that parton distribution functions 
defined for different n possess the same infrared structure, and thus the 
same evolution behavior, though different ultraviolet structure 0. After the 
derivation, we bring n back to the light cone, and the gauge invariance of 
a parton distribution function is restored ||. That is, the vector n appears 
only at the intermediate stage of the derivation, and acts as an auxiliary tool. 

The key step of the CS formalism is to obtain the derivative p + dF/dp + . 
Because of the scale invariance of F in the vector n as indicated by the gluon 
propagator, —iN^ u (l)/l 2 , with 

N = gT ; + n 2 —-- 2 , 2 

n ■ I (n ■ iy 

F must depend on p + via the ratio (p-n) 2 /n 2 . Hence, we have the chain rule 
relating p + d/dp + to d/dn [|], [ETJ: 



d „ n 2 d „ 

= v a — F , (3) 

dp + v ■ n dn a 
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with v a = 5 a+ a dimensionless vector along p. The operator d/dn a applies 
to a gluon propagator, giving @] 

d -N^ = l —(FN av + FN"") . (4) 



dn a n ■ I 

The loop momentum (l u ) in the above expression contracts with a vertex 
the differentiated gluon attaches, which is then replaced by a special vertex 
11 



n 2 v a 

v a = . 5 

v ■ nn ■ I 

This special vertex can be easily read off from the combination of Eqs. (0) 
and (§. 

Employing the Ward identities for the contraction of (l v ) ||, the sum 
of diagrams with different differentiated gluons reduces to a new diagram, in 



which the special vertex moves to the outer end of a parton line j|, |10[. We 
obtain the formula 

p + ^F(x,k T ,p + ) = 2F(x,k T ,p + ) , (6) 

shown in Fig. 1(a), where F denotes the new diagram, and the square rep- 
resents the special vertex. The coefficient 2 comes from the equality of the 
new functions with the special vertex on either of the two parton lines. 

The leading regions of the loop momentum flowing through the special 
vertex are soft and hard, since the vector n does not lie on the light cone, and 
collinear enhancements are suppressed. If following the standard resumma- 
tion procedure, we should factorize out subdiagrams containing the special 
vertex order by order in the leading regions (see Figs. 1(b) and 1(c)). This 
procedure will lead to a new unified evolution equation in the next section. 
To reproduce the CCFM equation, however, the factorization must be per- 
formed in a different way described by Fig. 2(a), where the two jet functions 
J group all-order virtual corrections, and the lowest-order real gluon emis- 
sion between them is soft. This factorization picture is motivated by the 
conventional derivation in [[3]. 

We first resum the double logarithms contained in J by considering its 
derivative 

P + ^J(Pt,P + ) = J(k T ,P + ) 
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= [Kj(k T //j,, a s (/j)) + Gj(p + /fi, a s (fj))] J(k T ,p + ) , (7) 

which is similar to Eq. (|6]). In the second expression subdiagrams containing 
the special vertex have been factorized into the functions Kj and Gj in the 
leading soft and hard regions, respectively. At lowest order, Kj comes from 
Fig. 2(b), where the replacement of the gluon line by an eikonal line will 
be explained below, and Gj comes from Fig. 2(c). The second diagram in 
Fig. 2(c) serves as a soft subtraction, which guarantees a hard momentum 
flow in Gj. We have set the infrared cutoff of Kj to k?, as indicated by its 
argument, in order to regularize its soft pole. This cutoff is necessary here 
due to the lack of the corresponding real gluon emission diagram. At this 
step, an ambiguity in the summation of single logarithms is introduced. 

We evaluate the contribution from Fig. 2(b). The color factor is extracted 
from the relation f a bcfbdc = —N c 5 a d, where the indices a,b, . . . are referred to 
Fig. 1(b), and N c = 3 is the number of colors. The one- loop Kj is written as 

r d 4 l yPY^ x 

Kj = iN ^iw? N ^^m- SKj > < 8) 

with 5Kj an additive counterterm. The triple-gluon vertex is 

r ^A = g^(xp + l) x + g vX (xp-2iy + g x ^(l-2xp) y , (9) 
« xp + (g^v x + g' /X v^-2g x W) , (10) 

where Eq. ( |10"D is the soft loop momentum (/ — > 0) limit of Eq. ([|). The factor 
l/(—2xp ■ I) is the soft approximation of the gluon propagator 1/(1 — xp) 2 . 

It has been shown that the terms v x and v M in Eq. (|10"D result in contribu- 
tions suppressed by a power 1/ s, s = (p + q) 2 , compared to the contribution 
from the last term v v [||, |6| . Absorbing the metric tensor g Am into F, Eq. (^) 
is simplified into 

The factor v u /(v ■ I) is represented by an eikonal line in Fig. 2(b), with the 
Feynman rules v v for an eikonal vertex and l/(v-l) for an eikonal propagator. 

Since ultraviolet poles exist in Fig. 2(b) and in the second diagram of 
Fig. 2(c), the sum Kj + Gj is independent of //, i.e., renormalization-group 
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(RG) invariant. The standard RG analysis leads to 



^T Kj = ~ lJ = ~^ Gj ' (12) 



whose solution is written as 

fp + dp 

Kj(hr/fJi,oi s {fji)) + Gj(p + /fi,a s (fi)) = -jj(a s ) / — . (13) 



The anomalous dimension jj is found to be 



7j = V-^&Kj = a s , (14) 

with the coupling constant a s = N c a s /n. To the accuracy of leading (double) 
logarithms 0, we have neglected the initial conditions Kj(l,a s (kT)) and 
Gj(l, a s (p + )) of the RG evolution, and the running of a s , which contribute 
to next-to-leading (single) logarithmic summation. In the derivation of the 
new unified evolution equation these factors can be included. 
Substituting Eq. ( |T3"D into Eq. (0), we solve for 

J(k T ,Q) = A 1 /\k T ,Q)J^ , (15) 

with the double-logarithm exponential 

r Q dp + f 'P + dp 



A(k T , Q) = exp 



-2a 



k T p^ Jk T /i 



(16) 



where the upper bound of p + has been chosen as Q. The initial condition 
can be regarded as a tree-level gluon propagator, which will be eikonalized 
in the evaluation of the lowest-order soft real gluon emission below. 
We split the above exponential into 

A{k T , Q) = A s (zq, Q)A NS (z, q, k T ) , (17) 

with z = x/£ and q = It/(1—z), where ^ is the momentum fraction entering J 
from the bottom, and It is the transverse loop momentum carried by the soft 
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real gluon. The so-called "Sudakov" exponential As and the "non-Sudakov" 
exponential A^s are given by 



A NS (z,q, k T ) 



exp 
exp 
exp 
exp 



-2cL 



-a. 



-2a, 



-a* 



Q dp + fP + dp 



zq Jk T \i 

Q 2 dp 2 ri-kr/p dz' 

(zq) 2 p 2 JO 1 — 2 

z i dp + fP + dp' 

k T p + Jk T ft 

k T/i dz' [ k2 T dp 



(18) 



Z 1 J{z'q) 2 p 2 



(19) 



where the variable changes p = (1 — z')p and p + = p for As, and p = p and 
p + = z'q for Atvs have been adopted to obtain the second expressions. 
Using Eq. (|l^) for the two jet functions, F from Fig. 2(a) is written as 



F(x,k T ,p + ) = iN c g 2 J ^LN uP (l)^ T 2m6(l 2 )A(k Tj Q) 



(2tt) 4 ~ '"" w vl 
x6(Q - zq)F(x + l + /p + , |k T + l T \,P* 



(20) 



where the tree-level gluon propagator on the right-hand side has been 
eikonalized to give the factor 1/v ■ I, and that on the left-hand side has been 
absorbed into F. The extra 9 function, requiring Q > zq, renders the Su- 
dakov exponential As meaningful. This relation comes from the angular 
ordering of radiative gluons, Q/(xp + ) > It/[{£ — x)p + }. Compared with the 
transverse momentum ordering for the DGLAP equation, Q (It) has been 
divided by the longitudinal momentum xp + ((£ — x)p + ). Therefore, the 
inserted scale zq reflects the specific kinematic ordering for the CCFM equa- 
tion. Those radiative gluons, which do not obey angular ordering, contribute 
to the non-Sudakov exponential. 

Adopting the variable change £ = x + l + /p + and performing the integra- 
tion over l~, we obtain 



F(x, kT,p^ 



a 



d 2 l 



2n 2 (£ — x)p~* 



vr [n+l 2 T + 2n-(i-x) 2 p+ 2 ] 2 ^ kT ' Q) 



x9(Q-zq)F(Z,\k T + l T \, P +) 



(21) 
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where n = (n + , n~ , 0) has been chosen for convenience. The above expression 
is then substituted into Eq. (§) to find a solution of F. Integrating Eq. (pip 
from p + = to Q, and applying the variable changes £ = x/z and It = 
(1 — z)q, we have 



1 /-7^ 

F{x,k T ,Q) = F^ + a s [ dz f -^e{Q-zq)As{Q,zq)A NS {z,q,h 

Jx J TTq 



X- 



z(l - z 



-F(x/z, \k T + (1 - z)q|,g) 



(22) 



where the initial condition F^ ' corresponds to the lower bound of p + . To 
workout thep + integration, F(x/z, |k^+lr|,p + ) in Eq. (pl|) has been approx- 
imated by F(x/z, |pt + 1t|j <?)• Note that the n dependence has disappeared, 
and the CCFM kernel is gauge invariant. 
Equation (|22"D can be reexpressed as 



F{x,k T ,Q) = 

with the function 
P = 



'dz / ^6(Q - zq)A s (zq,Q)P(z,q,k T ) 
xF{x/z,\k T + {l-z)q\,q) , (23) 



1 v + A iYS (^, g, fcr)- - 2 + z(l - ^) 
1 — ^)4. 2 



(24) 



which is close to the splitting function 



P, 



mi 



—!-— + --2 + z(l-z) 
l-z)+ z 



(25) 



Equation ( p3|) is the CCFM equation ||. To arrive at Eq. (|24|), we have 
employed the identity l/[z(l — z)} = 1/(1 — 2) + 1/z, and put in by hand 
the last term —2 + z(l — z) . This term, finite at z — > and at z — > 1, can 
not be obtained in either due to a soft approximation. Note that Eq. ([22]) 
contains an infrared singularity from z — > 1. Hence, we have replaced the 
factor 1/(1 — 2) by the plus distribution 1/(1 — z) + in Eq. (|23|), and dropped 
the soft pole corresponding to z — > 1. On the other hand, only the non- 
Sudakov form factor A^s in front of 1/z is retained. Because Ajv S vanishes 
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when the upper bound zq of p + approaches zero, as indicated by Eq. fllTf), it 
smears the z — > pole in P ||. 

As shown above, the double-logarithm exponential A and the function P 
are infrared divergent individually. To have a well-defined evolution equation, 
the infrared regularizations for A and P must be implemented. The lower 
bound of the variable z is x, and never goes down to zero. Therefore, the 
smearing factor A^s and the splitting of the double-logarithm exponential 
in Eq. flT7|) are not essential. The new unified evolution equation derived in 
the next section will take into account these observations. 



3. New Unified Equation 

In this section we shall derive a new unified evolution equation following 
the standard procedure of the CS formalism. We start with Eq. (||) for the 
unintegrated gluon distribution function F(x, kx,p + )- We Fourier transform 
Eq. into b space with b being the conjugate variable of k?- In the leading 
soft and hard regions of the loop momentum flowing through the special 
vertex, F can be factorized into the convolution of subdiagrams containing 
the special vertex with the original distribution function F. To lowest order, 
the soft contribution from the subdiagrams in Fig. 1(b) is written as F s = 
F sv + F sr with 

Fsv = iN c g 2 J -^N^(l)^F(x,b,p + ) , (26) 

F sr = tN c g" J ^N uP (l)^2nz5(iy^ b 

xF(x + l + /p + ,b,p + ) , (27) 

corresponding to virtual and real gluon emissions, respectively. Equation (p6|) 
is the same as the integral in Eq. (|TT|). The factor exp(zlr • b) in Eq. (f2"7j) 
comes from Fourier transformation of the real gluon contribution. 
We reexpress the function F in the integrand of F sr as 

F(x + l + /p + , b, p+) = 0((1 - x)p + - l + )F(x, b, p + ) 

+ [F(x + l + /p + ,b,p + ) -0((1 -x)p + -l + )F(x,b,p + )} . (28) 
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The contribution from the first term on the right-hand side of Eq. ( p8| ) is 
combined with F sv , leading to K(bfi,a s (fi))F(x,b,p + ) with the function 



K 



iN c g 2 
1 



dH 



■N vP (l) 



v f3 v u 



X 



, + 27ri5(l 2 )6((l - x)p + - l + )e ilT h 



5K 



(29) 



where 8K is an additive counterterm. The contribution from the second term 
in Eq. ( ^8|) is denoted as F' s (x,b,p + ), 

dH 



F' 



iN c g z 



2\ c U T b 



(2tt) 4 "'^ v ^ vl 
x [F(x + l + /p + , 6, p + )-6((l- x)p + - l + )F(x, 6, p + )} . (30) 



That is, we have separated the real gluon contribution into two pieces. The 
piece with an infrared pole, along with the virtual gluon contribution, goes 
into the function K. Compared with Kj in Eq. ([TT]) , K is infrared finite due 
to the inclusion of part of the real correction. The other piece F' s , being also 
infrared finite, will lead to a splitting function below. 

Fig. 1(c) gives G{k + / fi,a s (fi))F(x,b,p + ) with the function 



G = iN c g 2 



d 4 l AT vP 
(2^^ (/) ^ 



(/ - 2k + vY v L 



(I - k+v) 2 v ■ I 



-5G 



(31) 



where k + = xp + is the parton momentum and SG an additive counterterm. 
To obtain the above expression, Eq. (Q) for the triple-gluon vertex has been 
inserted, which generates the first term in the brackets. At intermediate x, 
G is characterized by a large scale k + , whose logarithms \nk + imply that 
collinear divergences are present in F, and that the Sudakov resummation 
of double logarithms from the overlap of collinear and soft enhancements 
is necessary. Through this function G, the dependence on a large scale is 
introduced into the new evolution equation as stated before. 

The functions K and G contain ultraviolet divergences individually as 
indicated by the renormalization scale fi. However, their ultraviolet diver- 
gences, both from the virtual contribution F sv , cancel each other. This is 
obvious from the counterterms 



6K 



-5G 



I + iln^-^, . 
e 2 2 ' 



(32) 
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with the Euler constant 7^. Therefore, the sum K + G is RG invariant, the 
same as Kj + Gj in Eq. (]T3|) . The results of K and G from Eqs. (|29|) and 
(§T]) are 



K(bn,a s {ij)) = -a s (ji) 
G(k + /[i,a s (fi)) = -a s (jj) 



K (2(l-x)p + ub)+\n b -^ + lE 



' 2A;+ 

in h In f 



(33) 



with the gauge factor v = 
The RG solution of K 



v ■ n 



l /\n 2 



G is given by 



K(b/i, a s (fi)) + G(k + /fi, a s (fi)) = 

K(l, a s (k + )) + a s (k + )) - s(b, k + ) 



with 



K(l,a s (k + )) = a s (k 

G(l,a s {k + )) = -a s {k + )\n2u 
r k+ dft 
n/b ft 



ln(l -x) + \n(p + b) + \n2is 



s(b, k + ) 



d 

1k{ois{^)) + /3(g)— K{1, a s (ji)) 



(34) 

(35) 
(36) 

(37) 



7if = fid5K/dfi being the anomalous dimension of K and (3(g) the beta 
function. To obtain Eq. (|35D, we have kept only the first two terms in the 
power series of the Bessel function K (x) ~ — ln(x/2) — j E , which is a rea- 
sonable approximation at small b. We shall show that the large b region is 
suppressed by the Sudakov form factor from double-logarithm resummation. 
Even if there is no Sudakov suppression, the intrinsic b dependence of the 
gluon distribution function F(x,b,p + ), usually parametrized as a Gaussian 
form |n| , still provides suppression. 

The one-loop expression of 7^ in Eq. ( |37D gives leading (double) logarith- 
mic summation, and the two-loop expression gives next-to-leading (single) 
logarithmic summation. Another term proportional to (3dK/dg and the ini- 
tial conditions K(l, a s ) and G(l, a s ) also give next-to-leading contributions. 
Hence, the Sudakov resummation can be performed up to the accuracy of 
single logarithms, and the running of the coupling constant must be taken 
into account in our formalism. Note that it is equivalent to choose Ck + 
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as the upper bound of Jjl in Eq. fl3~T|), where C a constant of order unity. 
The ambiguity from this arbitrary constant C is cancelled by the ambiguity 
from the initial condition G(l / C, a s ) at the single-logarithm level. The terms 
involving the gauge factor In 2v in the initial conditions of K and G cancel. 

For F' s , we employ the variable change £ = x + l + /p + , and perform the 
integration straightforwardly, obtaining 

F' s = 2a s {k + )p + vb J dftfi(2(f- x)p + vb 



a s (k 



F^,b,p + )-F(x,b,p A 

x 



i- x /i 



(38) 



For a similar reason, we have taken the small b limit of the Bessel function 
K\. At this stage, the gauge dependence from v in the evolution kernel 
disappears completely. The argument of a s in the above expression has 
been chosen as k + according to Eq. (|3"1D, since F' s is a part of the function 
K. Further applying the variable change £ = x/z and using the identity 
l/[z(l - z)\ = 1/(1 - z) + l/z, Eq. (M) is rewritten as 



F> = a s (k" 



1 J f 1 1 
dz(- + 

x \Z 1 — Z 



F(x/z, b, p + ) — F(x, b, p^ 



(39) 



Combining Eqs. ( |34"D and (|39|) , the evolution equation becomes 



p 



d 
dp^ 



F(x, b,p~* 



s(b, k + ) - a s (k + ) \n(k + b)} F(x, b,p + ) 



+2a s (k + ) [ dzP gg (z)F(x/z,b,p^ 

Jx 



(40) 



with P„„ defined by Eq. 



|). The term a s ln(l — x) in Eq. (|35| ) has been 
incorporated with 1/(1— z) in Eq. (^) to form the plus distribution 1/(1— z) + 
appearing in a DGLAP splitting function. The term — / dzF(x,b,p + )/ z in 
Eq. (^) has been combined with ln(p + b) in Eq. (p5|) , leading to ln(fc + 6) in 
Eq. ([5]. The terms finite at z — > and at z — > 1 in P gg are put in by 
hand as in derivation of the CCFM equation ||. In the derivation below we 
shall ignore the term a s (k + ) ln(k + b) in Eq. fl40|), which is less important than 
s(b, k + ). The integral of s gains a large contribution from the low end of the 
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integration variable p, ~ 1/6, since the running coupling constant a s (l/b) is 
much greater than a s (k + ) in the large 6 region. 

The Sudakov form factor A can be extracted by assuming 

F(x,b,p+) = A(x,b,p + ,Q )F^ s \x,b,p + ) , (41) 

with the exponential 

A(x,b,p + ,Q ) = exp 

where Q is an arbitrary low energy scale, satisfying xQ > 1/6. As xQ < 
1/6, implying the vanishing of collinear enhancements, the Sudakov form 
factor should be set to unity. Equation (|42|) indicates that the Sudakov form 
factor exhibits strong suppression at large 6, and that the small 6 expansion 
in the derivation of Eqs. ( |35|) and (|38|) makes sense. 

Substituting Eq. (pH] ) into Eq. (§3), we have the differential equation of 
± ) 

p+ -f«(,,6, p+) = Afe ;; p+ ; Qo) / i ".(w^pi . (43) 

The above equation is rewritten as 

F^(x, 6, Q) = F^(x, 6, Q ) + f dz f ^ ^^j^ F(x/z, 6, „) , (44) 

A Jq2 A(x,b,p,Q ) 

where F^ s \x,b,Qo) = F(x,b,Q ) is the initial condition of F at the low 
scale p + = Qo- We have chosen the upper bound of p + as the measured kine- 
matic variable Q. Finally, using Eq. (|4l|) , the unintegrated gluon distribution 
function in 6 space is given by 

F{x,b,Q) = A(x,b,Q,Q )F(x,b,Q ) 

rl rQ 2 d/J 2 _ 

+ dz —a s (xp)A(x,b,Q,p)P gg (z)F(x/z,b,p) , 

(45) 

Equation fl45T) is the new unified evolution equation. 



fp dp, . 
-2 / — s(b,xp) 
JQo p, 



(42) 
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4. Features of the New Equation 



In this section we shall explore features of the new unified evolution equa- 
tion. We first demonstrate that the new equation approaches the DGLAP 
equation at large x, The strong fey ordering of radiative gluons corresponds 
to small transverse separation b among these gluons. In the region with small 
b and large x, we have 

a s (xfj) « a s (ii) , (46) 
lim s(b,k + )=0, (47) 

6^1/fc+^0 

lim F(x,b,Q) =xG(x,Q) . (48) 

b— >0 

Equation ( |4"7| ) corresponds to the fact that the effect of Sudakov resummation 
is weak at small b as explained in Sec. 3. Equation fl48|) comes from the 
definition of the gluon density G(x, Q) as the integral of F(x, kr, Q) in Eq. (fj) 
over kx- 



f d 2 kT 

\imF(x,b,p + ) = / F(x,k T ,p + ) 



d 2 k 7 

71 

1 fdy 

xlY,(p,a\F+(y-,y T )F^(0)\p,*) 



p+ 



= xG(x,p + ) . (49) 
Differentiating Eq. (flop with respect to Q 2 under Eqs. (f46D-([48l), we obtain 



d r ' dz 

Q 2 ^G(x, Q) = a s (Q) / —P gg (z)G(x/z, Q) , (50) 

which is the DGLAP equation for G(x, Q). The initial condition F(x, b, Qq) 
does not depend on Q, and thus its derivative vanishes. 

We then show that the new evolution equation reduces to the BFKL 
equation in the small-x limit. For convenience, we shall adopt the one-loop 
expression j K = a s , which is the same as 7j for the CCFM equation, and 
ignore the running of a s below. Strong rapidity ordering corresponds to the 
assumption that the z integral in the second term of Eq. (|45|) is dominated 
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by the behavior of F{x/ z,b, fi) at small x/z, i.e., at 2; > x 0. Substituting 
F(x,b,Q) for F(x/z,b,Q), Eq. (]45D becomes 

F(x,b,Q) = A(x,b,Q,Q Q )F(x,b,Q ) 

ri dz fQ 2 da 2 
+a s / — / -^-A(x,b,Q,fj.)F(x,b,fj.) . (51) 

Jx Z JQq /i 

In the above expression the term —2 + z(l — z) in the splitting function P gg , 
which is finite — > 0, has been neglected. 
In the small x limit we have 



A(x,b,Q,Q )*l, (52) 

d Q 2 Q 2 

-x—A(x, b, Q, Qo) = a s In — ^ A(x, b, Q, Q ) w a s In — f . (53) 

ax C^o Vo 



Differentiating F with respect to x under the above relations, Eq. (^T|) leads 
to 



- x-^F(x,b,Q) = a s [ ( ^-F(x, b,fj) , 

ax JQl jjL A 

O 2 

« a s In ^ F(x, 6, Q) . (54) 

The derivative of the first term on the right-hand side of Eq. (^l|) has been 
neglected, since F^ is assumed to be "flat" |Tl[] . The derivatives of F(x, 6, //) 
and of A(x,b,Q, fi) in the integrand of the second term have been also ne- 
glected, because they give next-to-leading-order contributions. To derive the 
second expression, we have made the approximation F(x, b, //) ~ F(x,b,Q), 
which maintains the feature of the evolution in x. Equation ([54]) is the mod- 
ified BFKL equation with a Q dependence derived in ||. 

If we drop hard virtual gluon contributions, the Q dependence will not 
appear in the evolution kernel. Replacing Q by 1/6, which is the only scale, 
Eq. (|53]) reduces to 

- x— F(x, b) = -a s \n(Q 2 b 2 )F(x, b) . (55) 

(JjJb 

Since the evolution kernel does not depend on Q, the argument Q of F has 
been suppressed. It can be shown that the Fourier transformation of the 
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BFKL equation adopted in || [TT 
d „, . , r d 2 l 



x —F(x,k T ) = a s J -^[F(x,\k T + \ T \) - 6(Q - l T )F(x,k T )} , (56) 



coincides with Eq. (|55j). The conventional BFKL equation predicts that the 
small-x rises of the gluon distribution function and of structure functions 
involved in deep inelastic scattering are almost Q-independent as indicated 
by Eq. fl56f) . However, recent HERA data of the structure function F 2 (x, Q 2 ) 
[fLl^l exhibit a more sensitive Q dependence: the rise is rapider at larger Q. 
It is easy to observe that the evolution kernel in Eq. (|54l) enhances the rise 
of the gluon distribution function when Q is large. It has been shown that 
predictions from Eq. ( |54|) are well consistent with the data ||. 

It is known that the BFKL rise of the gluon distribution function and of 
the structure function F 2 violates the unitarity bound F 2 < const, x ln 2 (l/x). 
This is obvious from the solution of F to Eq. ([55]), which grows like a power of 
1/x. Without approximating F(x/z,b, /i) by F(x,b,fi), the term 1/(1 — z) + 
in P gg gives a destructive contribution to the evolution kernel because of 
F(x/z, b, /i) < F(x, b, /i). It implies that the assumption of rapidity ordering 
overestimates real gluon contributions, which are responsible for the BFKL 
rise. Without assuming strong rapidity ordering, Eq. (f55| ) becomes 

F{x,b,Q) = A{x,b,Q,Q )F{x,b,Q ) 

r 1 dz fQ 2 da 2 . . . „ ^ , , . 
+c *s — o — o, Q, H)F{x, b, ii) 
Jx z JQl a z 



'Qo W 

Q 2 

Z JQl fl* 



+ ^ -dz ^dj^ HxAQ ^ )F {x/zA»)-F{xA») 



(57) 

which contains a destructive term compared with Eq. (|5T|). Differentiating 
Eq. (|57D with respect to x and using Eq. (|53D , we derive 

-x-^F(x,b) = -a s ln(Qlb 2 )F{x,b) 



+ [a s HQ 2 b 2 )} 2 I" ^ g^jhZM , (58) 
J x z 1 — z 



dx 



where the variable Q has been replaced by 1/6 as in Eq. (155]). The correction 
term is of higher order, whose additional power a s ln(QQ& 2 ) comes from the 
derivative of A with respect to x. 
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Equation ([58]) is similar to the modified BFKL equation with unitarity 
proposed in 0. The behavior of F(x/z, b) at z — > x provides strong suppres- 
sion to real gluon contributions, such that the power-law rise turns into a 
logarithmic rise, and unitarity is restored j7|. The reasoning is reviewed be- 
low. Inserting a guess solution F oc x~ x into Eq. fl58D , A being a parameter, 
we obtain 

X = S + S 2 f 1 -^^, (59) 

Jx Z 1 — Z 

with S = — a s ln(Qo& 2 ). A solution of A, < A < S, exists for a x which is 
not too small, since the correction term, i.e., the second term on the right- 
hand side of Eq. fl59|), is negative and not large. That is, F increases as a 
power of x, consistent with the results from the conventional BFKL equation. 
While the correction term divere; GS clS X — > 0, and no solution of A is allowed, 
implying that F can not maintain the power-law rise at extremely small x. 
We then substitute another guess F\ x ^ oc ln(l/x) with a milder rise into 
Eq. (j58|). In this case the correction term, increasing as ln(l/x) roughly, 
cancels the first term SF oc ln(l/x). The derivative term —xdF/dx oc 1 
is negligible. Therefore, Eq. ([58]) holds approximately. At last, we assume 
F\ x _> oc const, as a test. It is easy to find that the first term becomes 
dominant, and the correction term and the derivative term vanish, i.e., no 
const. exists. These simple investigations indicate that F should increase 
as ln(l/x) at most, instead of as a power of 1/x, when x approaches zero. 
For a rigorous numerical solution to the modified BFKL equation, which 
exhibits a logarithmic rise at x — * explicitly, refer to [0. In conclusion, 
the new evolution equation, predicting a rapid power-like rise of the gluon 
distribution function F for small x and a milder logarithmic rise for x — > 0, 
satisfies the requirement of unitarity. 

We emphasize that predictions from the CCFM equation, based on strong 
angular ordering of radiative gluons, do not satisfy the unitarity bound. A 
careful investigation shows that the BFKL rise embedded in the CCFM equa- 



tion is attributed to the same origin as in the new evolution equation |6], [11 

a s C— I — 2 Q{Q - zq)A(k T , Q)F{x, |k T + (1 - z)q\,q) , (60) 
Jx z J irq z 

where As(zq, Q) and Ans(z, q, &t) have combined into A(kx,Q)- The cor- 
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rection to the power-law rise also arises from the term 1/(1 — z) + : 

- [ X dz f d 2 q F(x/z,\k T + (1 - z)q\,q) 

a ° ~ —0{Q ~ zq)A s {zq,Q) r . (61) 

Jx z J 7i q — + 

However, the Sudakov exponential A$(zq,Q) in Eq. (|6lD suppresses the re- 
gion with z — > x — > 0, such that the correction to the assumption of ra- 
pidity ordering in the CCFM equation is not important enough to turn the 
power-law rise into a logarithmic rise. This correction is not suppressed by 
A(x, b, Q, fx) at small x in the new equation. 

Compared with the CCFM equation, the new evolution equation bears 
other similarities and differences. They both contain the Sudakov expo- 
nentials, but the forms are different. The upper bounds of the evolution 
variables are the same (Q), but the lower bounds are zq in the CCFM equa- 
tion, and [L in the new equation. The longitudinal variable \i controls the 
magnitude of the Sudakov form factor by changing the partition of collinear 
enhancements between A(x, b, Q, //) and F(x/z, b, fi). While zq from angular 
ordering changes the partition of collinear enhancements between Ag(zq, Q) 
and A NS (z, q, There is no separation of A NS from As in the new equa- 
tion. As explained in Sec. 2, this smearing factor A^s is not necessary, since 
the z — > pole does not exist. 

5. Conclusion 

In this paper we have derived a new unified evolution equation for the 
gluon distribution function using the CS resummation technique, which is 
appropriate for both the large and small x. The features of this equation 
are summarized below. The infrared cancellation between real and virtual 
gluon emissions is explicit in both the Sudakov form factor and the splitting 
function. The Sudakov form factor can take into account next-to-leading 
logarithmic summation systematically, including the running of the coupling 
constant. Hard virtual gluon emissions introduces a desired Q dependence 
into the evolution kernel. The splitting function does not contain the non- 
Sudakov form factor compared to that in the CCFM equation. Predictions 
for the gluon distribution function satisfy the unitarity bound. Hence, the 
new unified evolution equation can be regarded as an improved version of 
the CCFM equation. 
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In the future studies we shall explore the properties of this new equation 
in details, and results will be published elsewhere. 

This work was supported by the National Science Council of R.O.C. under 
Grant No. NSC87-2112-M-006-013. 
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Figure Captions 



FIG. 1. (a) The derivative p + dF/dp + in the axial gauge, (b) The lowest- 
order subdiagrams for K. (c) The lowest-order subdiagrams for G. 

FIG. 2. (a) The subdiagram containing the special vertex for the CCFM 
equation, (b) The lowest-order subdiagram for Kj. (c) The lowest-order 
subdiagrams for Gj. 
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